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For QCD at non-zero chemical potential jx, the Dirac eigenvalues are scattered in the complex 
plane. We define a notion of ordering for individual eigenvalues in this case and derive the dis- 
tributions of individual eigenvalues from random matrix theory (RMT). We distinguish two cases 
depending on the parameter a = jX^F^V , where V is the volume and F is the familiar low-energy 
constant of chiral perturbation theory. For small a, we use a Fredholm determinant expansion and 
observe that akeady the first few terms give an excellent approximation. For large a, all spectral 
correlations are rotationally invariant, and exact results can be derived. We compare the RMT 
predictions to lattice data and in both cases find excellent agreement in the topological sectors 
v = 0,1,2. 
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1. Introduction 

Studies of the properties of the Dirac operator spectrum in gauge theories, including QCD, 
have a long history. For example, the low-lying Dirac modes provide information about sponta- 
neous chiral symmetry breaking through the Banks-Casher relation. The Dirac operator spectrum 
is also a natural object to study in lattice QCD. In the deep infrared, QCD in the £-regime can be 
described by the chiral random matrix theory (RMT) introduced in Ref. [1]. One of the advantages 
of RMT is that many exact analytical results can be derived. These results contain the low-energy 
constant (LEC) £ of chiral perturbation theory (chPT), and in some cases also the LEC F. These 
LECs can then be determined by fitting lattice data to the RMT curves. 

The observables that are most natural to compute in RMT are the spectral density correlation 
functions. At zero chemical potential pL, all of them are known in RMT [1, 2]. On the other 
hand, one can consider individual Dirac eigenvalue distributions (lED). From the lattice QCD point 
of view, these are the most natural observables to measure directly. Certain quantities such as 
the average positions of the eigenvalues are more pronounced in lEDs and therefore require less 
statistics to be measured reliably. 

At jLt = 0, the Dirac operator is anti-Hermitian and has a purely imaginary spectrum. In this 
case, all lEDs are known analytically in RMT [3]. They have become a standard tool in lattice QCD 
to extract £ in sectors of fixed topology. At /x / 0, the Dirac operator is no longer anti-Hermitian, 
and its eigenvalues are scattered in the complex plane. Our work is based on an RMT for /x 7^ 
[4] which has an eigenvalue representation and for which all complex density correlations (both 
quenched and unquenched) have been calculated. The same results can be obtained from the RMT 
for /I 7^ introduced earlier by Stephanov [5] or from chPT in the £-regime [6, 7] and are universal 
in that sense. These results have been compared to data from quenched lattice simulations with 
staggered [8] and overlap [9] fermions at /x / 0. A virtue of /i 7^ is that it couples to F in leading 
order of chPT [10] so that a comparison with lattice data allows us to extract F [8]. 

For the lEDs at /x 7^ much less is known. One of the problems here is to define an ordering of 
complex eigenvalues. Previous work on the repulsion between complex eigenvalues in RMT [1 1] 
and on the lattice [12] was done in the bulk of the spectrum, where no link to chPT is apparent. In 
the present work, we are interested in lEDs for eigenvalues close to the origin since they provide 
information on topological properties and LECs. 

We first define the general notion of lEDs for complex eigenvalues, and then compute the first 
few lEDs approximately by truncating a so-called Fredholm determinant expansion to the first few 
terms. It was already observed in Ref. [13] for /x =0 that this is a very good approximation. For 
large values of the parameter a = jJL^F^V , where V is the volume, we are able to derive all lEDs in 
closed form [14]. We use these exact results as a consistency check of the Fredholm determinant 
expansion. Our results are then compared to the lattice data of Ref. [9], in which the overlap Dirac 
operator for /x 7^ was constructed. This operator has good chiral properties (it satisfies a Ginsparg- 
Wilson relation, has an exact lattice chiral symmetry and exact zero modes, and satisfies the index 
theorem), which is essential for the present work. The same operator is obtained if a chemical 
potential is introduced in the domain-wall fermion formalism in the limit of infinite extent of the 
fifth dimension [15]. Due to the sign problem at /i 7^ 0, the lattice analysis in Ref. [9] was restricted 
to the quenched case, and this restriction on the lattice data applies to this work as well. 
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2. Individual eigenvalue distributions for complex eigenvalues 

Consider an operator with a finite number N of complex eigenvalues, distributed according to 
a joint probability distribution ^{{z}) which is symmetric in all its arguments. We also assume a 
z — > — z symmetry and only consider the upper half-plane C+. The partition function is then given 
by Z = J^ rijLi d^Zj^{{z}), and the spectral density correlation functions are defined as 

Rk{zu:.,Zk) = 2 (N-k)l / '^ ^k+i---d zn^n{zi,..;Zn)- (2.1) 

Now consider any one-parameter family of mutually non- intersecting closed contours dJ[ri] which 
cover C+. For fixed T], dJ[ri] is the boundary of a set /[tj]. Let us parametrize the contour as 
dJ{z{'c)) with z(t) = x{z) + iy{z). We then define the ^-th eigenvalue distribution pk{J, t) as the 
probability that k — l eigenvalues are inside /, one is at the point z(t) on the boundary dJ, and the 
remaining N — k aiein the complement /, 

p'^'^ ") " I G) J\I/'^ ,_n l/^i ^(i^i) L.^.(.) • (2.2) 

Note that the eigenvalue ordering is induced by the entire contour family. The {t] , t} play the role 
of generalized polar coordinates. In the following, the argument / of p^ will be suppressed. It is 
possible [14] to express all pk{'v) through the densities Eq. (2.1). In particular, for the distribution 
of the first eigenvalue one obtains 

pi{z) = RM^))-JdhiR2{zuz{z)) + ^^ J dhi j dh2R3{zuZ2,z{r)) + ... (2.3) 

One can show that the integrated distributions /\(^) = J dzj{ri,z)pk{ri,z) are normalized as 
In^!^ PkWdri = 1, where [T7min,T]max] is the range of T], and j{ri,z) = \{dx,dy)/{dri,dz)\ is the 
Jacobian of the transformation from {x,y) to (tj , t) . The proof is similar to that for real eigenvalues 
[13]. We emphasize that in this framework the choice of the contour family becomes part of the 
definition of the quantities we measure (i.e., the pts)- Different contour families in general lead to 
different p^s. However, one relation always holds trivially, namely Lf=i Pk{z) = Ri (z)- 

3. RMT predictions 

The partition function of the matrix model we use [4] reads 

Zv = [dAdB exp{-A^tr(A4l" +SS^)} f] det ( . ^ ""^'^ iA+fiB\ ^^^^ 

Here, A and B are complex (N + v) xN matrices with no further symmetries, v > is the topologi- 
cal charge, the m/ are the masses of Nf flavors of dynamical quarks, and /t is the chemical potential 
in the matrix model. In the laige-N limit this model describes QCD in the e-regime. All density 
correlation functions of this model follow from the kernel ^a?(z,,z*) of bi-orthogonal polynomials 
with respect to the weight 
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Figure 1: Top: Microscopic density pi(^) (left) and distributions p\{^) (middle) and p2{^) (right) of the 
first and second eigenvalue for an elliptic family of dJ[rj] parametrized as (Re^/5)^ + (Im^)^ — 11^ onC+, 
all for V = and a — 0.174. Bottom: Counterparts for ^ = 0. Left: Density. Middle: Density (dashed/blue), 
exact pi{£,) (solid/green), and approximate pi{£,) (dotted/red). Right: Density (long dashes/blue), exact 
pi{^) + pii^) (short dashes/purple), exact p2{^) (solid/green), and approximate p2{^) (dotted/red). The 
deviations of the approximate pi and p2 from the exact results are an artifact of the truncation of the Fred- 
holm determinant expansion. 



where Ky (and ly below) are modified Bessel functions, according to 

k 
Rk{zi,...,Zk 



Hw^^^'^Hz^) det KN{zi,z*j)=: det ^^(z,-,z}) 



(3.3) 



'r,2, 



VF^H^) 



We rescale tlie parameters of the model such that the parameters a = 2Nfi [= Vt^pi^), r\f = 
Nnif (= VLrrif), and i^jt = Nzk (= y^Zk) stay finite in the large-A^ limit. The scaling of these param- 
eters in terms of the LECs of chPT is given in parentheses. In the quenched case, the microscopic 
kernel J^(z,-,z*) = YimN^^J^Nizi/N ,z*/N)/N is given by [4] 
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and the microscopic spectral density follows as Pi((§) = J(fs{^^^*)- 

Approximate computations for arbitrary a 

In order to obtain the distribution p\, we use Eq. (2.3) and substitute the densities from 
Eqs. (3.3) and (3.4). Similar formulas can be written down for other pk&. For practical pur- 
poses, we truncate the (so-called Fredholm determinant) expansion in Eq. (2.3) to the first three 
terms and the coiTcsponding expansion for p2 to the first two terms. The density Pi(<^) and the 
distributions p\{^) and P2{^) of the first and second eigenvalue are shown in Fig. 1, together with 
their counterparts in the /i = case. For /i = exact results are available [3] which facilitate a 
detailed comparison. As for /i = [13], we see that the expansion converges rapidly. Higher-order 
terms merely assure that pk{^ ) remains zero for large \t,\. 
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Figure 2: Left: First eigenvalue distribution as a function of the radius, for the Ginibre ensemble and for 
the V = 1 chiral ensemble in the a ^f °° limit from Eq. (3.6). The v = 1 sector was chosen since for the 
comparison we need exactly one eigenvalue at the origin in both cases. Middle; Spectral density Eq. (3.5) 
and distributions of the first eight eigenvalues from Eq. (3.6) (and similar for pk{r) with fe > 2) as well as 
their sum, all in the a ^ °° limit and for v = 0. Right: Same for v = 1. 

Exact results in the large- a limit 

In the a ^ oo limit, the problem becomes radially symmetric. For finite but large a, the 
symmetry is apparent close to the origin. The limiting microscopic spectral density expressed in 
the new variable t, = (^/2\/a reads 



limpi(^: 



m 



n 



-K.m'Mm. 



(3.5) 



In this limit, we can derive a closed expression for all eigenvalue distributions [14]. Because of the 
rotational symmetry we choose dJ to be a semi-circle in C+ of radius r = |(^ | and obtain for pi (r) 



Pi{r) 



1 3^f r^^+^-+^^v+i(r^; ^^.2 
i=0 



TlrdrLM 2^^+^£\{£+v)\ 



■^v+i (r^)/lV? (^') + ^v+2(r2)C {r') 



(3.6) 



Here, we have introduced the incomplete Bessel function ly (x) = Y,i=o{^/'2.)^"^^ /n\{n + v)\ for 
^ > 0, and zero otherwise. Our result is analogous to the result for the nearest-neighbor spacing 
distribution of the Ginibre ensemble [11, 12], which can be inteipreted as the distribution of the 
smallest nonzero eigenvalue if one eigenvalue is fixed at zero. Expressions for Pk{r) with ^ > 2 are 
also available [14, 16]. Fig. 2 shows several p^, summing up nicely to the spectral density Eq. (3.5). 

4. Lattice calculations 



The lattice part of our work is based on the data obtained in Ref. [9]. The overlap Dirac 
operator introduced there is 

Dov(M) = l + 75£(75£»w(iU)), (4.1) 

where e is the sign function of a non-Hermitian matrix and Div(/i) is the Wilson Dirac operator at 
jU / 0. This overlap operator was shown to satisfy a Ginsparg-Wilson relation and to have good 
chiral properties [9, 15]. Equation (4.1) reduces to the standard overlap operator [17] at /x = 0. 

From the computational standpoint, the most demanding part is the computation of the matrix 
sign function. For the present set of data, this was done exactly using the spectral definition of the 
sign function. The lattice size is only V = 4^, since high statistics are needed for a comparison with 
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Figure 3: Integrated distribution P\ (R) of the first eigenvalue for v = 0, 1,2 and ^a = 0.1 (left), jxa = 0.3 
(middle), and /za = 1.0 (right). The solid lines are the RMT results from Eq. (2.3), the the histograms are 
the lattice data of Ref. [9]. The bending-up of the RMT curves for large R is an artifact of truncating the 
expansion (2.3). For jj.a— 1 .0 we also show the exact RMT results in the a ^ °° limit from Eq. (3.6). 

RMT. The coupling in the Wilson action is j3 = 5.1 in order to stay in the £-regime (where RMT 
applies) for the first eigenvalue(s) [9]. The Wilson mass is mwa = —2 (a is the lattice spacing), and 
the quark mass is zero. Data were sampled in the topological sectors v = 0, 1, 2 for the values of 
jXa = 0.1, 0.2, 0.3, and 1.0, coiTcsponding to a = 0.174, 0.615, 1.42, and 4.51, respectively. The 
number of configurations varied from about 9000 for ^ua = 0.1 to about 3000 for iJ.a = 1.0. The 
parameters £ and F were determined by a fit to the spectral density from Eq. (3.4). For iJ.a = 1.0, 
the data showed rotational invariance up to |^ | = 0.7, and hence only the combination L/F could 
be determined by a fit to Eq. (3.5). Our comparisons use these values and are thus parameter-free. 

To compute pi from the lattice, we choose for our contours dJ[ri] concentric semicircles with 
radius R = rj for all values of /i. (Other choices are also possible.) The localized nature of the 
lEDs allows us to integrate over the phase (i.e., we compute Pi (/?) = Jq ddRpi{R, 6)) rather than 
to consider cuts as in Ref. [9]. This procedure results in a much better signal. As a consequence, 
we are able to obtain a better comparison in topological sectors v = 0, 1 and, for the first time, to 
successfully test the RMT predictions in the v = 2 sector. In Fig. 3 we show the comparison of 
RMT predictions and lattice results for jxa = 0.1, 0.3, 1.0 and v = 0, 1,2 (the fia = 0.2 case is 
similar [16]). We see that for smaller /i, the agreement is excellent, whereas there are deviations 
for /la = 1.0 and v = 1, 2. In these cases we are outside the £-regime of QCD so that RMT no 
longer applies. We emphasize that while the ascent of the distributions from zero was in principle 
already tested in Ref. [9] through the density, their descent represents a new, parameter-free test. 

The deviations of the theoretical curves from zero for large R are an artifact of the trunca- 
tion of the Fredholm determinant expansion. As jj. (or a) is increased, the convergence of the 
approximation becomes slower, i.e., more terms are needed. 

For jUa = 1 .0 we are almost in the radially symmetric regime. Thus we expect Pi (7?) computed 
approximately through Eq. (2.3) to be close to the exact result Eq. (3.6) in the a ^ oo limit. This 
expectation is confirmed in Fig. 3 (right). 

5. Conclusions 



In this work we have studied individual Dirac eigenvalue distributions in the e-regime of QCD 
at nonzero chemical potential. We provided a general framework for ordering complex eigenvalues. 
Our RMT computation for arbitrary a was based on the truncation of a Fredholm determinant 
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expansion. In the a ^ oo limit, we were able to derive all lEDs analytically in closed form. These 
predictions were then tested against lattice data based on the generalization of the overlap Dirac 
operator to nonzero chemical potential, in the topological sectors v = 0, 1,2. We found excellent 
agreement between RMT and lattice results for several values of /x in the domain of the applicability 
of RMT. The descent of the lEDs represents a new, parameter-free test of RMT predictions. The 
much improved signal (resulting from the integration over the phase) allowed us, for the first time, 
to successfully test the RMT predictions in the topological sector v = 2. 

Acknowledgments 

This work was supported by EPSRC grant EP/D031613/1 (GA & LS), by EU network EN- 
RAGE MRTN-CT-2004-005616 (GA) and by DFG grant FOR 465 (JB & TW). 

References 

[1] E. V. Shuryak and J. J. M. Verbaarschot, Nucl. Phys. A 560 (1993) 306 [hep-th/9212088]; 
J. J. M. Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531 [hep-th/9401059]. 

[2] G. Akemann, P H. Damgaard, U. Magnea and S. Nishigaki, Nucl. Phys. B 487 (1997) 721 
[hep-th/9609174]; 
P. H. Damgaard and S. M. Nishigaki, Nucl. Phys. B 518 (1998) 495 [hep-th/971 1023]. 

[3] T. Wilke, T Guhr and T. Wettig, Phys. Rev. D 57 (1998) 6486 [hep-th/9711057]; 

S. M. Nishigaki, P H. Damgaard and T. Wettig, Phys. Rev. D 58 (1998) 087704 [hep-th/9803007]; 
R H. Damgaard and S. M. Nishigaki, Phys. Rev. D 63 (2001) 045012 [hep-th/00061 1 1]. 

[4] J. C. Osborn, Phys. Rev. Lett. 93 (2004) 222001 [hep-th/0403 1 3 1 ] ; 

G. Akemann, J. C. Osborn, K. Splittorff and J. J. M. Verbaarschot, Nucl. Phys. B 712 (2005) 287 
[hep-th/0411030]. 

[5] M. A. Stephanov, Phys. Rev. Lett. 76 (1996) 4472 [hep-lat/9604003]. 

[6] K. Splittorff and J. J. M. Verbaarschot, Nucl. Phys. B 683 (2004) 467 [hep-th/03 10271]. 

[7] F. Basile and G. Akemann, arXiv:0710.0376 [hep-th]. 

[8] G. Akemann and T Wettig, Phys. Rev. Lett. 92 (2004) 102002 [Erratum-ibid. 96 (2006) 029902] 
[hep-lat/0308003]; 
J. C. Osborn and T. Wettig, PoS LAT2005 (2006) 200 [hep-lat/0510115]. 

[9] J. Bloch and T. Wettig, Phys. Rev. Lett. 97 (2006) 012003 [hep-lat/0604020] . 

[10] D. Toublan and J. J. M. Verbaarschot, Int. I. Mod. Phys. B 15 (2001) 1404 [hep-th/0001 1 10]. 

[11] R. Grobe, R Haake and H.-I. Sommers, Phys. Rev. Lett. 61 (1988) 1899. 

[12] H. Markum, R. Pulhrsch and T. Wettig, Phys. Rev. Lett. 83 (1999) 484 [hep-lat/9906020] . 

[13] G. Akemann and P H. Damgaard, Phys. Lett. B 583 (2004) 199 [hep-th/03 11 171]. 

[14] G. Akemann and L. Shifrin (2007), to be pubhshed. 

[15] J. Bloch and T. Wettig, arXiv:0709.4630 [hep-lat], to appear in Phys. Rev. D. 

[16] G. Akemann, I. C. R. Bloch, L. Shifrin and T. Wettig, arXiv:07 10.2865 [hep-lat]. 

[17] R. Narayanan and H. Neuberger, Nucl. Phys. B 443 (1995) 305 [hep-th/9411108]; 
H. Neuberger, Phys. Lett. B 417 (1998) 141 [hep-lat/9707022]. 



